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ABSTRACT 

We present a method for defining higher-order moments of a spin-2 field on the sky using the trans- 
formation properties of these statistics under rotation and parity. For the three-point function of the 
cosmic shear we show that the eight logically possible combinations of the shear in three points can be 
divided into two classes, four combinations are even under parity transformations and four are odd. We 
compute the expected value of the even parity ones in the non-linear regime using the halo model and 
conclude that on small scales of the four combinations there is one that is expected to carry most of the 
signal for triangles close to isosceles. On the other hand, for collapsed triangles all four combinations are 
expected to have roughly the same level of signal although some of the combinations are negative and 
others positive. We estimate that a survey of a few square degrees area is enough to detect this signal 
above the noise at arc minute scales. 

Subject headings: cosmology: theory - gravitational lensing - large scale structure of universe- cosmic 
microwave background 



1. INTRODUCTION 

Gravitational lensing is quickly becoming a major tool 
for doing cosmology. Mass determinations of cluster of 
galaxies are now routine, there are several detections of 
weak lensing by the large scale structure of the universe 
and measurements of galaxy-galaxy lensing. Detections of 
weak lensing (Bacon et al. 2000; Kaiser et al. 2000; Van 
Waerbeke et al. 2000; Wittman et al. 2000) are partic- 
ularly important as they will provide constraints on the 
matter budget of the universe and an independent deter- 
mination of the power spectrum of the dark matter fluc- 
tuations over a wide range of scales that can be then com- 
pared with other measurements such as those from galaxy 
clustering or the anisotropics in the Cosmic Microwave 
Background (CMB). 

As a result of non-linear gravitational evolution, the pro- 
jected mass density field on the sky (k), which is responsi- 
ble for the deflections that lead to the measured shear, is 
expected to be highly non-Gaussian even if the initial seeds 
of density perturbations were perfectly Gaussian. The 
aim of observations that attempt to make projected mass 
maps is not only to measure the projected mass power 
spectrum, but also higher-order moments which contain 
additional information. It has been stressed that valu- 
able constraints on cosmological parameters are expected 
to come from a joint measurement of the variance and 
higher-order moments such as the skewness of the weak 
lensing convergence field because the variance is strongly 
dependent on both the amplitude of the mass power spec- 
trum and the density of the Universe, while the skewness 
(properly normalized) is essentially a measure of the latter 
(Bernardeau, Van Waerbeke & Mellier 1997; Jain & Seljak 
1997; Schneider et al. 1998; Van Waerbeke, Bernardeau & 
Melher 1999; Van Waerbeke et al. 2000). 

The non-linear growth of structure is very important on 
the length scales relevant for weak lensing, so estimates 
for the level of non-Gaussianity in k come from two differ- 



ent techniques: N-body simulations (Couchman, Barber 
& Thomas 1999; Jain, Seljak & White 2000; White & 
Hu 2000) and semi-analytic models (Jain & Seljak 1997; 
Schneider et al. 1998; Van Waerbeke et al 2001). N-body 
simulations have been used to create mock k maps from 
which higher-order moments have been measured. Semi- 
analytic techniques such as the halo model (see Cooray 
& Sheth 2002 for a recent review) or proposals about 
the behavior of higher-order correlations in the non-linear 
regime (Scoccimarro & Frieman 1999) have also been used 
to make weak lensing predictions (Hui 1999; Cooray, Hu 
& Mhalda-Escude 2000; Cooray & Hu 2001). 

So far most of these theoretical approaches have dealt 
with k; however the quantity that is directly measured is 
the shear 7^ . The crucial difference between k and 7 is that 
the shear is a spin-2 field on the sky while k is just a scalar. 
Thus the shear has two components at each position on 
the sky so when designing any statistics for the shear one 
has to make sure that the statistic does not depend on 
some arbitrary choice of coordinate system. Moreover the 
transformation from shear to k cannot be done exactly 
if one only has a shear map over a finite region on the 
sky and only at the position of the background galaxies. 
To circumvent this issue, measures of the shear at differ- 
ent points on the sky can be combined to form statistics 
that are invariant under rotation and can be expressed as 
weighted averages of k; the best example being the aper- 
ture mass (Kaiser et al. 1994, Schneider, van Waerbeke, 
Jain, & Kruse 1998). We will show that the amplitude of 
the three-point function varies with triangle configuration 
and can be both positive or negative. As a result not all 
combinations of the three-point function are optimal from 
a signal to noise perspective. 

A natural way to look for non-Gaussianity is then to look 

^What is actually measured is the reduced shear 7/(1 — k), thus 
on sufficiently small scales this may change theoretical predictions. 
The spin properties of the three-point function we will discuss apply 
both to the shear and the reduced shear 
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at statistics of the aperture mass. In this paper we will 
focus on another approach, directly defining higher-order 
statistics in terms of the shear field that are independent 
of the coordinate system. Any statistic of the shear can 
ultimately be written as a linear combination of the statis- 
tics we present in this paper. We will show however that 
not all the possible higher-order moments are expected to 
have the same level of cosmological signal and, moreover, 
the sign of these statistics can be positive or negative de- 
pending on the configuration of the points. The advantage 
of our approach is that it will allow us to isolate the con- 
figurations that have most cosmological signal and avoid 
suppression of signal to noise in the measurements or can- 
cellations that would result from an arbitrary linear com- 
bination. 

A first attempt has been made to define a three-point 
function for the shear field in Bernardeau et al (2002a). 
Their prescription corresponds to integrating over a par- 
ticular linear combination of the four different shear three- 
point functions that wc define in this work. They have ap- 
plied their statistic to a cosmic shear survey and reported 
a detection at approximately 5<t level (Bernardeau et al. 
2002b). Given these encouraging results it is worth to 
consider in more detail how to construct shear three-point 
functions and what is expected theoretically about their 
order of magnitude and sign depending on the particular 
configuration of the points. 

This paper will be written using weak lensing language, 
but identical issues arise if one wants to define higher-order 
moments of the CMB polarization field. The analogue of 
the shear components are the Q and U Stokes parameters 
and the analogue of the projected mass density is usually 
called the E field. Even if the initial conditions arc Gaus- 
sian, higher-order moments of the CMB polarization can 
be generated by secondary effects such as lensing so the 
statistics presented here will be equally applicable for the 
CMB. 

2. DEFINING HIGHER-ORDER MOMENTS FOR SPIN-2 
FIELDS 

In this section wc will show how to define higher-order 
moments of the weak lensing shear or CMB polarization 
fields in a way that is geometrically meaningful. The prob- 
lem is that the shear is a spin-2 field and thus at each point 
on the sky it has two components. Just as in the case of 
a vector field the values of these components depend on 
the choice of coordinate system. If at any given point one 
rotates the coordinate system used to define the shear com- 
ponents by an angle a (anticlockwise in our convention) 
the shear field has to be transformed as, 

j[ = cos 2a 7i -|- sin 2a 72 

72 = — sin 2a 71 -|- cos 2a 72 (1) 

Any meaningful statistic of the shear measured in a set 
of points has to depend only on the distances and relative 
orientation of these points and not on their absolute posi- 
tion on the sky or their orientation with respect to some 
fiducial origin. For example the two point function has to 
depend only on the distance between the two points and 
the three-point function only on the size and shape of the 
triangle formed by the three points. 

The way to overcome this problem of definition for the 



second moment is well known, one uses the separation vec- 
tor to define a "natural" coordinate system. The idea is 
to align the coordinate system so that one of the axis lies 
on the great circle joining the two points and uses that 
coordinate system to define the two components of the 
shear (Miralda-Escude 1991, Kaiser 1992, Kamionkowski, 
Kosowsky & Stebbins 1997). 

In general the way to define higher-order moments that 
are invariant under rotation is by contracting the measured 
shear at the three points with suitable combinations of the 
vectors that define the sides of the triangle, are invariant 
under translation and have the correct spin to compen- 
sate for the spin of the shear. Clearly this procedure is 
not unique as there are many such combinations. In this 
paper we will propose a simple, intuitive but at the same 
time geometrically meaningful way to define higher-order 
moments. 

An N-point function is characterized by the set of points 
X = {8,} {i = I, ... N) where the shear or the CMB 
polarization is measured. We denote 2D vectors on the 
sky with boldface. We then define the "center of mass" of 
X, 

1 ^ 

and use o as the origin when defining the shear at each of 
the points in X. At every point we can define a component 
of the shear along the direction that separates o and 6i 
which we can call 7+ and a component which is measured 
in the coordinate system that is rotated by 45°, which 
we call 7x . Figure 1 illustrates how the 7_|_ and 7x are 
defined. This is totally analogous to the procedure for the 
two point function, but now the two points used in the 
definition of 7+ and 7x are o and Oi instead of the two 
points in the two point function. The key to our procedure 
is that by defining an origin based only on the points in X 
we make sure that our statistic depends only on intrinsic 
properties of X and not on a fiducial origin. 

Another important property of spin-2 fields is that 
they can be decomposed into two scalar potentials, one 
that is even under parity (E) and one that is odd (B) 
(Kamionkowski, Kosowsky & Stebbins 1997, Zaldarriaga 
& Seljak 1997, Crittenden et al. 2001, Schneider et al. 
2002). In fact for finite sky coverage or for maps with 
holes, there is a third family of modes, the ambiguous 
modes for which one does not have enough information to 
decide whether the contribution is coming from E or from 
B (Lewis, Challinor, & Turok 2002, Bunn, Zaldarriaga, 
Tegmark, & de Oliveira- Costa 2002). Weak gravitational 
lensing only produces E modes {E is nothing but the pro- 
jected mass density k). 

The properties under parity transformation of 7+ and 
7x are different. As is clear from the definition, to obtain 
7x one has to rotate the coordinate system anticlockwise. 
This rotation changes direction when we do a parity trans- 
formation. This means that under parity 7^ — > 7+ but 
7'x ^ -7x- 

As a consequence of the difference in behavior of 7+ 
and 7x any estimator that contains an odd number of 7x 
is odd under parity. The difference in their parity be- 
havior will make some of the three-point functions van- 
ish. To illustrate what this means we can consider the 
case of the three-point function for equilateral triangles. 
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Fig. 1. — The point o is used to define the 7+ and 7x components 
of the shear at point 0. The two rods indicate the ellipticity at 
that would produce a positive 7_|_ and 7x 

In principle there are 8 different combinations of the two 
shear components at the three points, however only four 
of them: (7+7+7+), ( 7+7x7x ), ( 7x7+7x ), (7x7x7+) 
contain any signal from weak lensing because they are the 
only four that are even under parity. We show patterns 
that produce positive values for these correlations in Fig. 2. 
For a discussion on how the difference in behavior under 
parity of the three-point functions of general configura- 
tions constraints their possible values see Takada & Jain 
2002b. 

Being able to separate combinations that have signal 
from those that do not is very important not to dilute the 
signal one is trying to measure when combining or binning 
the measured statistics to obtain a detection. Moreover 
for equilateral triangles for example the four odd combi- 
nations: (7x7x7x), (7+7+7x>, (7x7+7+), (7+7x7+) 
may prove useful to monitor systematic problems in the 
data or to identify interesting physical effects such as clus- 
tering of the background sources, the effects of multiple 
scatterings, etc. (Bernardeau 1998, Jain, Seljak & White 
2000). 

Any connected N-point function of the shear can be 
written as a linear combination of the connected N-point 
functions of the two potential E and B. It is worth noting 
that in the case of connected N-point functions with odd 
number of legs there will be some configurations (such as 
those where the lenghts of all sides are equal) that can be 
separated into a set that only depends on the higher-order 
correlations of E and a set that only depends on the higher 
moments of B (assuming E and B are independent). We 
can argue this only on the basis of their bahavir under par- 
ity. The product of an odd number of i?'s is even under 
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Fig. 2. — The top row (from left to right) shows four shear patterns 
that would produce a positive value for (7+7+7+), (7x7x7+)i 
(7x7+7x {7+7x7x )• The bottom row shows four shear patterns 
that would produce a positive value for (7x7x7x ), (7+7+7x ), 
(7+7x7+ >. (7x7+7+ > 



parity while a product of an odd number of i?'s is odd, 
thus the even parity N-point functions of the shear receive 
contributions only from E and the odd ones only from B. 
For an N-point function with an even number of legs this 
is not true because the product of an even number of Es, 
or Bs is even under parity. The clearest example of this 
is the two-point correlation function where ( 7+7+ ) and 
(7x7x ) receive contributions from both {EE) and {BB) 
while (7x7+ ) is zero if there is no -B — -B cross-correlation. 

As we mentioned above the way to make a three-point 
function (or an N-point function for that matter) that is 
scalar is to contract the shear at the three points with some 
combination of the vectors that form the sides that trans- 
forms appropriately under rotations to cancel the spin of 
shear. That is we need to construct spin— 2 combinations 
of these vectors. Our proposed scheme is easy to under- 
stand. For each member of X we define 6 = 6 — and we 
construct two spin— 2 quantities, 

p+ = {el ~ 91,29 Jy)/e'' 

= (cos2(/), sin 20) 

P^ = {-2eJy,el-el)/9^ 

= (-sin 20, cos 2^). (3) 

The statistics we proposed are obtained by contracting the 
above quantities with the shear three-point function. For 
example, 

(7+7+7+) = ^r^f (7^,7;..7;.3), (4) 

where the index fi runs over the two components of the 
shear. The other three-point functions that we define are 
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obtained by replacing some of the P+ by Px • Finally we 
note that the vector di is nothing but 6i ~ [{Oi — 62) + 
{di — 03)]/3, i.e. basically the sum of the vectors that 
define the sides of the triangle that cross at 9\ . The same 
is true for the other vertices. 

3. WORKED EXAMPLES 

Our objective in this section is to get some intuition into 
how these higher-order correlations behave. To keep things 
simple we will just focus on the throe-point function. To 
get some idea of how these functions behave we will work 
in the context of the halo model, i.e. the dark matter is 
assumed to be distributed in a collection of halos of dif- 
ferent mass. Although analytic approximations and fits to 
numerical simulations exist for the profile of these halos 
and their mass function, this modeling is clearly simplistic 
to model the shear field as in reality halos are not spheri- 
cal. This particularly is bound to affect the configuration 
dependence of the three-point function. Thus definite pre- 
dictions will need more detailed modeling and comparison 
with direct measurements using numerical simulations. 

Our aim in this paper is more modest, we only want 
to gain some insight into how these different three-point 
functions behave, if they are positive or negative for exam- 
ple. We will start by considering the simple case of lensing 
by a singular isothermal sphere. We do this because the 
calciilation of the three-point function can be done ana- 
lytically, and it provides a useful check on the numerical 
code used to do the same calculations in the halo model, 
presented in 3.2. 

3.1. The Singular Isothermal Sphere 

In this section we will calculate the three-point function 

that results from a single halo with a power-law density 
run. For definiteness we will write down formulas for a 
singular isothermal sphere (SIS), but other power laws can 
be calculated in analogous way. For a SIS the density 
p{r) depends on distance from the center r as p(r) oc r~^; 
therefore, both the projected mass density n and the shear 
7 scale with projected separation as rj^. 

The shear pattern around a spherical halo is tangential 
centered at the origin of the halo which we will call u. 
When defining the shear components at a point 0, we will 
need to rotate the shear elements so as to define them 
relative to the vector 6 — o. We will call this angle of 
rotation a. We can write 



cos(2a) 
|6>-u| 



7xW = 



sin(2Q!) 

i6»-ur 



(5) 



The cosine and sine of 2a can be calculated in terms of 

{e-u)-{e- o) 



C0s(q;) 
sin (a) 



\e-u\ \0-o\ 

z ■ [(0 - u) X (6) - o) 

\e-u\ \e-o\ 



(6) 
(7) 



where z is the unit vector perpendicular to the plane of the 
sky. The three-point function ( is obtained by integrating 
over the position of the center of the halo u, 



, J s f ,2 cos(2q;i 



cos(2q;i) cos(2q;2) cos(2q:3) 



|6»2-u| \e.-n\ 



(8) 



c+""ai,/2,;3) 



cos(2ai) sin(2a2) sin(2a3) 

\di — u| |02 — u| |03 — u| 



,(9) 



where h^hyh sue the lengths of the triangle sides. If ~ 
|6>2 - e>i|2, q = 1^3 - 02p, I3 = \0i - ^sP- There are 
two permutations of the second equation which give the 
remaining three-point correlators, and C7 

To calculate Eqs. (8-9), we proceed as follows. First, 
for simplicity, we take the origin of coordinates to coincide 
with the center of mass, so o vanishes. We then write the 
u dependence inside the cosine and sine in Eqs. (6-7) in 
terms of di = Oi — u. This can be done simply using that 



U di 



(10) 



whereas the magnitude of u can be conveniently written 
as 

= i(d? + dl + 4) -^{11 + il + il). (11) 

Similarly we have = (2^| + 2?f - ;|)/9 and cycHc 

permutations. In this way, after appropriate translation 
in u the integrals in Eqs. (8-9) are of the form. 



J(l/l,J/2,J^3) 



{V?Y^ [(ll-u)2]-2[(l2_u)2]^ 



(12) 



which can be evaluated by using dimensional regular- 
ization techniques (see Scoccimarro 1997 and references 
therein) in terms of Apell's hypergeometric function of two 
variables, F4, with the series expansion: 



=0 j=o 



{c)i{d)j 



(13) 



where {a)i = r{a + i)/r{a). In our case, the arguments of 
F4 have a special symmetry that allows us to write them 
in terms of regular hypergeometric functions, 

F4a,j + j' -a- l,j,j';x{x- l),y{y- 1)] = 
F{a, 7 7' - a - 1, 7, x) X F{a, 7 7' - a - 1, 7', y), 

(14) 

which can be easily evaluated on the computer using 
MATHEMATICA. Figure 3 shows the results for different 
triangle configurations. Since the SIS is scale-free, the 
overall size of the triangle scales the results by a factor, 
so we take h = 1. The top panel shows the four (^'s 
for I2 = Zi as a function of the angle </> between Ii and 
I2. We see that for most of the triangles (^ ^ 0.37r), the 
signal is dominated by C^^"*", with a maximum close to 
equilateral triangles (0 = 27r/3)^. On the other hand, the 
three contributions involving 7x are generally smaller in 
magnitude. 

This can be understood geometrically from Fig. 4. Let 

us consider the halo center to be inside the triangle, which 
minimizes the distances to the vertices and thus maximizes 
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Fig. 3. — The shear three-point functions ^7(^11 '2, '3) for the case 
of isothermal sphere halos, as a function of the angle <j) between li 
and I2 for I2 = h (top panel) and I2 = 2l\ (bottom panel). Different 
line styles correspond to the four shear correlators, as labeled in the 
top panel. 

the signal. Since cos(2a) > for \a\ < ir/i, as long as 
the internal angles of the triangle are smaller than tt/2, 
is positive. As (j) ^ tt, a2 ^ 0, whereas |ai| = jaal 
cannot be larger than tt/2; thus C.^^^ remains positive 
as (j) n. On the other hand, as (/) — > 0, ai = ^3 — > 0, 
whereas \a2 \ — > tt if u is off center; that explains why C-^"*"*" 
becomes negative as (f> ^ 0. 

For the three other three-point functions involving 7x , 
things are more subtle. Since sin(2a) changes sign at 
a = 0, 7x can be positive or negative depending on the 
location of u relative to the bisector of the vertex. As 
a result moving u inside the triangle leads to cancella- 
tions, and thus smaller amplitude for Ct's- The bottom 
panel in Fig. 3 shows results for configurations in which 
I2 = 2li. A similar pattern is seen, where for isosceles 
triangles ~ O.Gtt the C.^"*"*" is positive and maximum, 
whereas the remaining three contributions are smaller, be- 
coming comparable only for collapsed triangles which can 
avoid cancellations. 

The physical interpretation of the different amplitudes 
is sketched in Fig. 5. In the top (bottom) row we show 
patterns that would produce positive (negative) values for 
C:J"'"+, C7 C^"*"^; and C^^^ from left to right, respec- 
tively. It is clear from the figure that the top examples 
correspond to patters that would be produced by overden- 
sities and the bottom patterns would be produced by un- 
derdensities. We indicate with a square the region where 
the overdensity (underdensity) should be located to pro- 
duce such pattern. 

3.2. Superposition of NFW profiles 

In this section we will calculate the three-point func- 
tions using the halo model (Peacock & Smith 2000, Sel- 




FlG. 4. — Definition of triangle variables. We characterize the 
shear generated by a halo with center at u by measuring its com- 
ponents with respect to the center of mass of the triangle, denoted 
by o. 




Fig. 5. — Examples of shear patterns generated at the vertices 
of equilateral triangles by overdensities (top panel) and underdensi- 
ties (bottom panel) located at the position denoted by square sym- 
bols. From left to right, contributions to Ct"^"*"! C-y ^"''i C7 
(^j^^ (vertices labeled anticlockwise starting from bottom-left ver- 
tex). Note that on the left- most plots, the density perturbation is 
located at the center of the triangle. 

jak 2000, Ma & Fry 2000, Scoccimarro et al. 2001). For 
simplicity we will restrict ourselves to the one-halo term 
which dominates on the small scales where the three-point 
function of the shear is easiest to measure in observations. 
For examples of calculations of weak gravitational lensing 
higher-order moments in the context of the halo model see 
Cooray & Hu 2001 and Takada & Jain 2002a. Measure- 
ments of higher-order moments of the convergence field in 
numerical simulations are given in Jain, Seljak & White 
(2000), White & Hu (2000); Van Waerbeke et al (2001) 
also present results for aperture mass statistics which is 
directly related to the cosmic shear. 

Under our assumptions the averaged shear can be writ- 
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ten as an integral over radial distance (dx = c dt/ a), mass 
(M) and angular location (u) of halos, 



(7(^1,^2,^3)= / dxd\{x) / dM 



dn 



7(u,0i,M,x) 7(u,02,M,x) 7(u,03,M,x), (15) 

where a is the expansion factor of the universe, dA is 
the comoving angular diameter distance, dn/dM is the 
mass function of halos. We have introduced the notation 
j{u,6, M,x) to indicate the shear produced at position 
by a halo at radial distance x angular position u. 
Note that for convenience we are using the same symbol 
regardless of whether 7+ or 7x is involved. 

To obtain an estimate for the three-point function we 
will evaluate Eq. (15) assuming that the background sources 
used to measure the shear are all at redshift Zg = I and 
that the cosmological model is the so called LCDM model 
( = 0.3, Ha = 0.7, h = 0.7, erg = 0.9 and n = 1). We 
will assume that the mass function of halos is that given 
by Sheth & Tormen 1999, Jenkins et al. 2001 and that 
dark matter halos have an NFW profile (Navarro, Frenk 
& White 1997). In particular, we use 



p{r) = 



Ps 



r/rs{l + r/rsf 



(16) 



where r is measured in comoving coordinates and is 
related to the virial radius of the halo by the concentra- 
tion parameter c, = Tmr/c. The mass of the halo is 
given by M = A■Kpsrl^^{\xl{l + c) - c/(l + c))/c?. The 
virial radius is calculated using that M = 47r/3r3^^poA(z) 
with po the mean density of the universe today and A(z) 
the overdensity of collapse as a function of redshift (ie. 
A(z = 0) « 340 for LCDM). For the concentration we 
take c{M,z) = 9/(1 + z){M/M^)-°-^^ (Bullock et al. 
2001), where M* is the mass contained in a sphere of 
radius i?* at which the variance of the density is one 
(i?* « 3.14 h-^Mpc for LCDM). 

The shear produced by an NFW profile at the origin is 
calculated as follows. The shear is expressed in terms of 
second derivatives of the gravitational potential -0, 71 = 
l/2(ipxx — 4'yy) 72 = ijjxy The gravitational potential 
satisfies 



k{u) 



2k(u) 
1 



dz 



pi^u^d\{x)+z^) (17) 



where 1/Ecrit = ATrGdAix)dA{x - Xs)/c^dA{Xs), X gives 
the radial position of the halo and Xs that of the back- 
ground sources^. Equation (17) is easily integrated to ob- 
tain ij} because k is only a function of u. Once the shear is 
obtained for a halo at the origin, we obtain 7(u, 0,M, %) 
with a coordinate transformation. A similar evaluation 
for the SIS to compare with the results obtained by the 
method of the previous section gives a useful check to our 
numerical integration code. 

Figure 6 shows the results of our calculation for some 
specific triangles. In the top panel h ^ I2 = 2' and in the 
bottom panel the case li — 2' and I2 ~ 4'. The similari- 
ties with the results for the SIS are striking, thus we can 

^In this equation c is the speed of light not the concentration 
parameter of dark matter halos 
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Fig. 6. — Even components of the three-point function of the shear 
calculated using the 1-halo term for LCDM. The top (bottom) panel 
shown the case = /2 = 2' (^i = 2', I2 = 4'). 



understand the dependence of each of the curves with in 
exactly the same way. 

The expected level of C7 on arcminute scales is C,^ ~ 
10~^. We estimate the number of triangles {Nt) needed 
to detect the three-point function in a particular configu- 
ration above the noise produced by the intrinsic ellipticity 
of the background galaxies as follows. If the typical ellip- 
ticity of the background galaxies is ~ 0.3 then the typical 
noise added to each components of the intrinsic shear is 
(7 ~ 0.3/\/2 ~ 0.2. The expected noise in the three-point 
function is then o'^/^/Nt- This implies that one needs 
roughly Nt ^ 6.5 x 10^ triangles of a particular shape to 
estimate the three-point function. In a survey covering a 
solid angle with a mean desity of galaxies n there are 
roughly Nt = nVl (TrR^n)^ /6 triangles with sides of scale 
R. Current surveys have about 20 galaxies/arcmin^ imply- 
ing that in a few square degrees there are enough triangles 
with sides of order an arcminute to detect this signal. 

4. SUMMARY AND DISCUSSION 

In this paper we have introduced a new way of defining 
higher-order correlation functions of a spin- 2 field such as 
the weak lensing shear or the CMB polarization by using 
the "center of mass" of the configuration as the origin from 
which the components of the shear are defined. In princi- 
ple for an A''-point function there are 2^ different statistics 
of the shear. We have shown that these statistics can be 
divided according to their behavior under parity transfor- 
mation and that one does not expect a cosmological signal 
in the odd ones for some configurations, such as for equi- 
lateral triangles. 

In order to gain intuition about the behavior of these 
statistics we calculated the four even three-point functions 
under some simple assumptions. We calculated analyti- 
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cally what would be expected to be produced by an ensem- 
ble of singular isothermal spheres. We showed that C^^^ 
is positive and expected to carry the bulk of the signal for 
triangles that are not too elongated. For elongated tri- 
angles we showed that all four three-point functions have 
similar values but some of them are positive and others 
are negative. If two of the points are very close to each 
another C'^''""'' and C7 """^ carry most of the signal. 

We estimated the three-point functions in the context 
of the halo model using the contributions from the one- 
halo term, which should be a reasonable approximation at 
small angular scales. We showed that the configuration 
dependence in this case is almost identical to that found 
for the SIS. We estimated that in order to detect a signal 
above the noise (due to the intrinsic ellipticity of galaxies) 
at scales of order of one arcminute a survey of a few square 
degrees is necessary. 

Clearly both the fact that we restricted ourselves to 
the one-halo term and that we assumed the halos to be 
spherical will affect the behavior of the three-point func- 
tion. A detailed study using numerical simulations will 
be needed to improve upon the calculation presented here 
(Benabed, Scoccimarro & Zaldarriaga. in preparation). 
The potential rewards of detecting a non-Gaussian signal 
in the shear maps are enourmous and the tantalizing de- 
tections reported so far (Bernardeau, Mellier & Van Waer- 
beke 2002b) make this a very exciting time to study these 
issues in detail. 

When this paper was under completion, a similar pro- 
posal for calculating the shear three-point function was 
put forward by Schneider & Lombardi (2002). 
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ABSTRACT 

We present a method for defining higher-order moments of a spin-2 field on the sky using the trans- 
formation properties of these statistics under rotation and parity. For the three-point function of the 
lensing siiear we show that of the eight logically possible combinations of the shear in three points only 
four are expected to have cosmological signal. We compute the expected value of these statistics in the 
non-linear regime using the halo model and conclude that on small scales of the four combinations there 
is one that is expected to carry most of the signal for triangles close to isosceles. On the other hand, for 
collapsed triangles all four combinations are expected to have roughly the same level of signal although 
some of the combinations are negative and others positive. We estimate that roughly 10^ triangles are 
needed to detect this signal above the noise. 

Subject headings: cosmology: theory - gravitational lensing - large scale structure of universe- cosmic 
microwave background 



1. INTRODUCTION 

Gravitational lensing is quickly becoming a major tool 
for doing cosmology. Mass determinations of cluster of 
galaxies are now routine, there are several detections of 
weak lensing by the large scale structure of the universe 
and measurements of galaxy-galaxy lensing. Detections 
of weak lensing (Bacon et al. 2000; Kaiser et al. 2000; 
Van Waerbeke et al. 2000; Wittman et al. 2000) are 
particularly important as they will provide constraints on 
the matter budget of the universe and an independent 
determination of the power spectrum of the dark mat- 
ter fluctuations over a wide range of scales that can be 
then compared with other measures such as those coming 
from galaxy clustering or the anisotropics in the Cosmic 
Microwave Background (CMB). 

As a result of non-linear gravitational evolution, the pro- 
jected mass density field on the sky (/t), which is responsi- 
ble for the deflections that lead to the measured shear, is 
expected to be highly non-Gaussian even if the initial seeds 
of density perturbations were perfectly Gaussian. The 
aim of observations that attempt to make projected mass 
maps is not only to measure the projected mass power 
spectrum, but also higher-order moments which contain 
additional information. It has been stressed that valu- 
able constraints on cosmological parameters are expected 
to come from a joint measurement of the variance and 
higher-order moments such as the skewness of the weak 
lensing convergence field because the variance is strongly 
dependent on both the amplitude of the mass power spec- 
trum and the density of the Universe, while the skewness 
(properly normaHzed) is essentially a measure of the latter 
(Bernardeau, Van Waerbeke & Mellier 1997; Jain & Seljak 
1997; Schneider et al. 1998; Van Waerbeke, Bernardeau & 
MeUier 1999; Van Waerbeke et al. 2000). 

The non-Hnear growth of structure is very important on 
the length scales relevant for weak lensing, so estimates for 
the level of non-Gaussianity in the k come from two dif- 
ferent techniques: N-body simulations (Couchman, Bar- 



ber & Thomas 1999; Jain, Seljak & White 2000; White & 
Hu 2000) and semi-analytic models (Jain & Seljak 1997; 
Schneider et al. 1998; Van Waerbeke et al 2001). N-body 
simulations have been used to create mock k maps from 
which higher order moments have been measured. Semi- 
analytic techniques such as the halo model (see Cooray 
& Shctli 2002 for a recent review) or proposals about 
the behavior of higher-order correlations in the non-linear 
regime (Scoccimarro & Frieman 1999) have also been used 
to make weak lensing predictions (Hui 1999; Cooray, Hu 
& Miralda-Escude 2000; Cooray & Hu 2001). 

So far most of these theoretical approaches have dealt 
with k; however the quantity that is directly measured is 
the shear 7. The crucial difference between k and 7 is that 
the shear is a spin-2 field on the sky while k is just a scalar. 
Thus the shear has two components at each position on the 
sky so when designing any statistics for the shear one has 
to make sure that the statistic does not depend on some 
arbitrary choice of coordinate system. Moreover the trans- 
formation from shear to k cannot be done exactly if one 
only has a shear map over a finite region on the sky and 
only at the position of the background galaxies. To circum- 
vent this issue, measures of the shear at different points on 
the sky can be combined to form statistics that are scalars 
under rotation and can be expressed as weighted averages 
of k; the best example being the aperture mass (Kaiser et 
al. 1994, Schneider, van Waerbeke, Jain, & Kruse 1998). 

A natural way to look for non-Gaussianity is then to look 
at statistics of the aperture mass. In this paper we will 
focus on another approach, directly defining higher-order 
statistics in terms of the shear field that are independent 
of the coordinate system. Any statistic of the shear can 
ultimately be written as a Hnear combination of the statis- 
tics we present in this paper. We will show however that 
not all the possible higher-order moments are expected to 
have a cosmological signal and, moreover, the sign of these 
statistics can be positive or negative depending on the con- 
figuration of the points. The advantage of our approach 



1 



2 



is that it will allow us to isolate the configurations that 
have cosmological signal and avoid suppression of signal 
to noise in the measurements or cancellations that would 
result from an arbitrary linear combination. 

A first attempt has been made to define a three-point 
function for the shear field in Bernardeau et al (2002a). 
Their prescription corresponds to integrating over a par- 
ticular linear combination of the four different shear three- 
point functions that we define in this work. They have ap- 
plied their statistic to a cosmic shear survey and reported 
a detection at approximately 2(t level (Bernardeau et al. 
2002b). Given these encouraging results it is worth to 
consider in more detail how to construct shear three-point 
functions and what is expected theoretically about their 
order of magnitude and sign depending on the particular 
configuration of the points. 

This paper will be written using weak lensing language, 
but identical issues arise if one wants to define higher-order 
moments of the CMB polarization field. The analogue of 
the shear components are the Q and U Stokes parameters 
and the analogue of the projected mass density is usually 
called the E field. Even if the initial conditions are Gaus- 
sian, higher-order moments of the CMB polarization can 
be generated by secondary effects such as lensing so the 
statistics presented here will be equally appHcable for the 
CMB. 

2. DEFINING HIGHER-ORDER MOMENTS FOR SPIN-2 
FIELDS 

In this section we will show how to define higher order 
moments of the weak lensing shear or CMB polarization 
fields in a way that is geometrically meaningful. The prob- 
lem is that the shear is a spin-2 field and thus at each point 
on the sky it has two components. Just as in the case of 
a vector field the values of these components depend on 
the choice of coordinate system. If at any given point one 
rotates the coordinate system used to define the shear com- 
ponents by an angle a (anticlockwise in our convention) 
the shear field has to be transformed as, 

7^ = cos 2a 7i -I- sin 2a 72 

72 = — sin 2a 71 -I- cos 2a 72 (1) 

Any meaningful statistic of the shear measured in a set 
of points has to depend only on the distances and relative 
orientation of these points and not on their absolute posi- 
tion on the sky or their orientation with respect to some 
fiducial origin. For example the two point function has to 
depend only on the distance between the two points and 
the three point function only on the size and shape of the 
triangle formed by the three points. 

The way to overcome this problem of definition for the 
second moment is well known, one uses the separation vec- 
tor to define a "natural" coordinate system. The idea is to 
align the coordinate system so that one of the axis lies on 
the great circle joining the two points and uses that coor- 
dinate system to define the two components of the shear 
(Kaiser 1992, Kamionkowski, Kosowsky & Stebbins 1997). 

In general the way to define higher order moments that 
are scalar under rotation is by contracting the measured 
shear at the three points with suitable combinations of the 
vectors that define the sides of the triangle, are invariant 
under translation and have the correct spin to compen- 
sate for the spin of the shear. Clearly this procedure is 
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Fig. 1. — The point o is used to define the 7+ and 7x components 
of the shear at point 6. The two rods indicate the ellipticity at 
that would produce a positive 7-1- and 7x 

not unique as there are many such combinations. In this 
paper we will propose a simple, intuitive but at the same 
time geometrically meaningful way to define higher-order 
moments. 

An N-point function is characterized by the set of points 
X = {Oi} (i = 1, ... N) where the shear or the CMB 
polarization is measured. We denote 2D vectors on the 
sky with boldface. We then define the "center of mass" of 
X, 

1 ^ 

and use o as the origin when defining the shear at each of 
the points in X. At every point we can define a component 
of the shear along the direction that separates o and 6i 
which we can call 7+ and a component which is measured 
in the coordinate system that is rotated by 45", which 
we call 7x. Figure 1 illustrates how the 7-1- and 7x are 
defined. This is totally analogous to the procedure for the 
two point function, but now the two points used in the 
definition of 7+ and 7x are o and 0, instead of the two 
points in the two point function. The key to our procedure 
is that by defining an origin based only on the points in X 
we make sure that our statistic depends only on intrinsic 
properties of X and not on a fiducial origin. 

Another important property of spin-2 fields is that 
they can be decomposed into two scalar potentials, one 
that is even under parity (E) and one that is odd (B) 
(Kamionkowski, Kosowsky & Stebbins 1997, Zaldarriaga 
& Seljak 1997). In fact for finite sky coverage or for maps 
with holes, there is a third family of modes, the ambiguous 
modes for which one does not have enough information to 
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Fig. 2. — The top row (from left to right) shows four shear patterns 
that would produce a positive value for ( 7_|_7_|.7_|_ ), (7x7x7+ )i 
(7x7+7x )) (7+7x7x )• The bottom row shows four shear patterns 
that would produce a positive value for (7x7x7x )i (7+7+7x }: 
(7+7x7+ >: (7x7+7+ > 

decide whether the contribution is coming from E or from 
B (Lewis, ChalUnor, & Turok 2002, Bunn, Zaldarriaga, 
Tegmark, & de OHveira-Costa 2002). Weak gravitational 
lensing only produces E modes {E is nothing other than 
the projected mass density k). 

The properties under parity transformation of 7+ and 
7x are different. As is clear from the definition, to obtain 
7x one has to rotate the coordinate system anticlockwise. 
This rotation changes direction when we do a parity trans- 
formation. This means that under parity 7^ 7+ but 
7'x -7x- 

As a consequence of the difference in behavior of 7+ 
and 7x any estimator that contains an odd number of 7x 
is odd under parity. In the case of weak lensing we do not 
expect any signal in these configurations. 

To illustrate what this means we can consider the case 
of the three-point function. In principle there are 8 differ- 
ent combinations of the two shear components at the three 
points, however only four of them: { 7_|_7_|_7_|_ ), (7_|_7x7x ), 
{7x7+7x ), (7x7x7+ ) contain any signal from weak lens- 
ing because they are the only four that are even under 
parity. We show patterns that produce positive values for 
these correlations in Fig. 2. 

Being able to separate combinations that have signal 
from those that do not is very important not to di- 
lute the signal one is trying to measure when combin- 
ing or binning the measured statistics to obtain a detec- 
tion. Moreover the four odd combinations: (7x7x7x ), 
{7+7+7x), (7x7+7+), (7+7x7+) may prove useful to 
monitor systematic problems in the data or to identify 
interesting physical effects such as clustering of the back- 
ground sources. 

It is also worth noting that in the case of connected N- 



point functions with odd number of legs one can always 
separate them into a set that only depends on the higher- 
order correlations of E and a set that only depends on 
the higher moments of B (assuming E and B are inde- 
pendent) because the product of an odd number of Es is 
even under parity while a product of an odd number of Bs 
is odd. Under this circumstance the even parity N-point 
function receive contributions only from E and the odd 
ones from B. For even number of legs this is not true be- 
cause the product of an even number of Es or Bs is even 
under parity. The clearest example being the correlation 
function where (7+7+ ) and (7x7x ) receive contributions 
from both {EE) and (BB) while (7x7+ ) is zero if there 
is no E — B cross-correlation. 

As we mentioned above the way to make a three point 
function (or an N-point function for that matter) that is 
scalar is to contract the shear at the three points with some 
combination of the vectors that form the sides that trans- 
forms appropriately under rotations to cancel the spin of 
shear. That is we need to construct spin —2 combinations 
of these vectors. Our proposed scheme is easy to under- 
stand. For each member of X we define 9 = 9 — o and we 
construct two spin —2 quantities, 

p. = {-2eJy,el-el)/0'. (3) 

The statistics we proposed are obtained by contracting the 
above quantities with the shear three point function. For 
example, 

(7+7+7+) = PfPfPf (7MaM.7M3), (4) 

where the index n runs over the two components of the 
shear. The other three point function that we defined 
are obtained by replacing some of the P+ by Px • Finally 
we note that the vector 9i is nothing other that 9i = 
[{9i — ©2) + (^1 — ^3)1/3, basically just the sum of the 
vectors that define the sides of the triangle that cross at 
9i. The same is true for the other vertices. 

3. WORKED EXAMPLES 

Our objective in this section is to get some intuition into 
how these higher-order correlations behave. To keep things 
simple we will just focus on the three-point function. To 
get some idea of how these functions behave we will work 
in the context of the halo model, i.e. the dark matter is 
assumed to be distributed in a collection of halos of dif- 
ferent mass. Although analytic approximations and fits to 
numerical simulations exist for the profile of these halos 
and their mass function, this modehng is clearly simplistic 
to model the shear field as in reality halos are not spheri- 
cal. This particularly is bound to affect the configuration 
dependence of the three-point function. Thus definite pre- 
dictions will need more detailed modeling and comparison 
with direct measurements using numerical simulations. 

Our aim in this paper is more modest, we only want to 
gain some insight into how these different three-point func- 
tions behave, if they are positive or negative for example. 
We will start by considering the simple case of lensing by a 
singular isothermal sphere. We do this because the calcu- 
lation of the three-point function can be done analytically. 
Then we will consider the results in the halo model. 
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3.1. The Singular Isothermal Sphere 

In this section we will calculate the three-point func- 
tion that results from halos with a power-law density run. 
For definiteness we will write down formulas for a singu- 
lar isothermal sphere (SIS), but other power laws can be 
calculated in analogous way. For a SIS the density p{r) de- 
pends on distance from the center r as p{r) oc r~^; there- 
fore, both the projected mass density k and the shear 7 
scale with projected separation as rj^. 

The shear pattern around a spherical halo is tangential 
centered at the origin of the halo which we will call u. 
When defining the shear components at a point 9, we will 
need to rotate the shear elements so as to define them 
relative to the vector 9 — o. We will call this angle of 
rotation a. We can write 



1+iO) 



cos(2q;) 



sin(2Q;) 



(5) 



The cosine and sine of 2a can be calculated in terms of 

(e - u) - (0 - o) 



COs(q!) 

sin(a) 



\9-u\ \9-o\ 
z-[{9-u)x{9- o) 



\9- 



\9- 



(6) 
(7) 



where z is the unit vector perpendicular to the plane of the 
sky. The three-point function C is obtained by integrating 
over the position of the center of the halo u. 



cos(2q;2) cos(2q;3) 



|6>2 -u| |6>3 -u| 



(8) 



(9) 



A-hxx/, , ,^_ f ,2„ cos(2ai) sin(2a2) sin(2a3) 
^ i 1^1 - u| \92 - u| \93 - u| 

where h,l2,h are the lengths of the triangle sides, Ij = 
|02 - I'l = \93 - 02p, ll = 1^1 - There are 

two permutations of the second equation which give the 
remaining three-point correlators, C^"^^ and C7 

To calculate Eqs. (8-9), we proceed as follows. First, 
for simplicity, we take the origin of coordinates to coincide 
with the center of mass, so o vanishes. We then write the 
u dependence inside the cosine and sine in Eqs. (6-7) in 
terms of di = 9i — u. This can be done simply using that 



n.9, = l{ 



+ \9i 



0. 



(10) 



whereas the magnitude of u can be conveniently written 
as 

n' = l{di + 4 + 4)-l{q + q + il). (ii) 

Similarly we have \9i\'^ = (2/2 -|- 2ll - lj)/9 and cyclic 
permutations. In this way, after appropriate translation 
in u the integrals in Eqs. (8-9) are of the form. 



(m2)-i[(1i-u)2]-2[(12-u)^ 



(12) 



which can be evaluated by using dimensional regular- 
ization techniques (see Scoccimarro 1997 and references 



therein) in terms of Apell's hypergeometric function of two 
variables, F4, with the series expansion: 



F4{a,b;c,d;x,y) = 

i=0 j=0 



x'yi {a)i+j{b) 



i+j 



(c)i(d)j 



(13) 



where (a)j = T{a + i)/T{a). In our case, the arguments of 
F4 have a special symmetry that allows us to write them 
in terms of regular hypergeometric functions, 

F4[a,j + j' - a- l,j,j';x{x - l),y{y - 1)] = 
F(a, 7-1-7' -a- l,'j,x) X ^(0,7-1-7' -a- l,'j',y), 

(14) 

which can be easily evaluated on the computer using 
MATHEMATICA. Figure 3 shows the results for different 
triangle configurations. Since the SIS is scale-free, the 
overall size of the triangle scales the results by a factor, 
so we take li = 1. The top panel shows the four (^j's 
for I2 = /i as a function of the angle (j> between li and 
I2. We see that for most of the triangles (0 ^ O.Stt), the 
signal is dominated by (^^^^, with a maximum close to 
equilateral triangles (0 = 27r/3)^. On the other hand, the 
three contributions involving 7x are generally smaller in 
magnitude. 

This can be understood geometrically from Fig. 4. Let 
us consider the halo center to be inside the triangle, which 
minimizes the distances to the vertices and thus maximizes 
the signal. Since cos(2a) > for \a\ < 7r/4, as long as 
the internal angles of the triangle are smaller than 7r/2, 
C^~^~^ is positive. As — > tt, 02 ^ 0, whereas |q;i| = jaal 
cannot be larger than 7r/2; thus C^"*"*" remains positive 
as ^ TT. On the other hand, as ^ 0, ai = as — )■ 0, 
whereas |q;2 | — > tt if u is off center; that explains why C^"*"*" 
becomes negative as ^ ^ 0. 

For the three other three-point functions involving 7x , 
things are more subtle. Since sin(2Q;) changes sign at 

= 0, 7x can be positive or negative depending on the 
location of u relative to the bisector of the vertex. As 
a result moving u inside the triangle leads to cancella- 
tions, and thus smaller amplitude for C-/'^- The bottom 
panel in Fig. 3 shows results for configurations in which 
I2 =2/1. A similar pattern is seen, where for isosceles 
triangles (f) « O.Gtt the C^"*"*" is positive and maximum, 
whereas the remaining three contributions are smaller, be- 
coming comparable only for collapsed triangles which can 
avoid cancellations. 

The physical interpretation of the different amplitudes 
is sketched in Fig. 5. In the top (bottom) row we show 
patterns that would produce positive (negative) values for 
C^^~^, C-^"*"^) and Cj^^ from left to right, respec- 
tively. It is clear from the figure that the top examples 
correspond to patters that would be produced by overden- 
sities and the bottom patterns would be produced by un- 
derdensities. We indicate with a square the region where 
the overdensity (underdensity) should be located to pro- 
duce such pattern. 

3.2. Superposition of NFW profiles 

^The divergence of C'y'^'^ as c/i — f tt is a peculiarity of the SIS 
when (3 
suggests. 



0. 



is regular in that limit despite what Fig. 3 
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Fig. 3. — The shear three-point functions CtC^Ii ^2, ^s) for the case 
of isothermal sphere halos, as a function of the angle </) between li 
and I2 for ^2 = h (top panel) and 1-2, = 2li (bottom panel). Different 
line styles correspond to the four shear correlators, as labeled in the 
top panel. 





Fig. 5. — Examples of shear patterns generated at the vertices 
of equilateral triangles by overdensities (top panel) and underdensi- 
ties (bottom panel) located at the position denoted by square sym- 
bols. From left to right, contributions to t^^"'""'", Cj^'^: Cj'^^ : 
C^^^ (vertices labeled anticlockwise starting from bottom-left ver- 
tex). Note that on the left-most plots, the density perturbation is 
located at the center of the triangle. 




Fig. 4. — Definition of triangle variables. We characterize the 
shear generated by a halo with center at u by measuring its com- 
ponents with respect to the center of mass of the triangle, denoted 
by o. 



In this section we will calculate the three-point func- 
tions using the halo model (Peacock & Smith 2000, Sel- 
jak 2000, Ma & Fry 2000, Scoccimarro et al. 2001). For 
simplicity we will restrict ourselves to the one-halo term 
which dominates on the small scales where the three-point 
function of the shear is easiest to measure in observations. 
For examples of calculations of weak gravitational lensing 
higher-order moments in the context of the halo model see 
Cooray & Hu 2001 and Takada & Jain 2002. Measure- 
ments of higher-order moments of the convergence field in 
numerical simulations are given in Jain, Seljak & White 
(2000), White & Hu (2000); Van Waerbeke et al (2001) 
also present results for aperture mass statistics which is 
directly related to the cosmic shear. 

Under our assumptions the averaged shear can be writ- 



ten as an integral over radial distance {dx = c dt/a), mass 
(M) and angular location (u) of halos, 

(7(^1, '2,^3) = I dxd\{x) I dM-^ I d'u 

7(u,0i,M,x) j(u,02,M,x) j(u,es,M,x), (15) 

where a is the expansion factor of the universe, dA is 
the comoving angular diameter distance, dn/dM is the 
mass function of halos. We have introduced the notation 
■y{u,6, M,x) to indicate the shear produced at position 
by a halo at radial distance x ^md angular position u. 
Note that for convenience we are using the same symbol 
regardless of whether 7+ or 7x is involved. 

To obtain an estimate for the three-point function we 
will evaluate Eq. (15) assuming that the background sources 
used to measure the shear are all at redshift Zg = 1 and 
that the cosmological model is the so called LCDM model 
( n„ = 0.3, = 0.7, h = 0.7, erg = 0.9 and n = 1). We 
will assume that the mass function of halos is that given 
by Sheth & Tormen 1999, Jenkins et al. 2001 and that 
dark matter halos have an NFW profile (Navarro, Frenk 
& White 1997). In particular, we use 



p{r) 



r/rs(l -I- r/r^y 



(16) 



where r is measured in comoving coordinates and is 
related to the virial radius of the halo by the concentra- 
tion parameter c, = r^ir/c. The mass of the halo is 
given by M = 47rpsr3 .^(ln(l + c) - c/(l -I- c))lc'. The 
virial radius is calculated using that M = 4ir/3rl^j,poA{z) 
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Fig. 6. — Even components of the three point function of the shear 
calculated using the 1-halo term for LCDM. The top (bottom) panel 
shown the case h = l2 = 2' (h — 2', h — 4'). 



with po the mean density of the universe today and A{z) 
the overdensity of collapse as a function of redshift (ie. 
A(z = 0) « 340 for LCDM). For the concentration we 
take c{M,z) = 9/(1 + z){M/M^)-°-^^ (Bullock et al. 
2001), where M* is the mass contained in a sphere of 
radius J?* at which the variance of the density is one 
(i?* « 3.14 h-^Mpc for LCDM). 

The shear produced by an NFW profile at the origin is 
calculated as follows. The shear is expressed in terms of 
second derivatives of the gravitational potential ip, = 
l/2{ipxx ~ i'vv) ci-iid 72 = ifxy The gravitational potential 
satisfies 



dz p{^u^d\{x)+z-^) (17) 



where l/Scrit = 47rG'd^(x)d.i(x - Xs) I c^dA^Xs), X gives 
the radial position of the halo and Xs that of the back- 
ground sources^. Equation (17) is easily integrated to ob- 
tain -0 because k is only a function of u. Once the shear is 
obtained for a halo at the origin, we obtain 7(u, 6, M, x) 
with a coordinate transformation. A similar evaluation 
for the SIS to compare with the results obtained by the 
method of the previous section gives a useful check to our 
numerical integration code. 

Figure 6 shows the results of our calculation for some 
specific triangles. In the top panel h = h = 2' and in the 
bottom panel the case li = 2' and 1-2 =4'. The similari- 
ties with the results for the SIS are striking, thus we can 
understand the dependence of each of the curves with 4> in 
exactly the same way. 

The expected level of on arcminute scales is ~ 

^In this equation c is the speed of light not the concentration 
parameter of dark matter halos 



10~^. We estimate the number of triangles (Nt) needed 
to detect the three point function in a particular configu- 
ration above the noise produced by the intrinsic ellipticity 
of the background galaxies as follows. If the typical value 
of one of the components of the intrinsic shear is a ~ 0.1, 
the expected noise in the three point function is cf^/^/Nt- 
This impHes that one needs roughly Nt ~ 10^ triangles 
of a particular shape to estimate the three-point function. 
Current surveys have about 20 galaxies / arcminute^ im- 
plying that in a few square degrees there are enough tri- 
angles to detect this signal even when divided into a few 
different triangle configurations. 

4. SUMMARY AND DISCUSSION 

In this paper we have introduced a new way of defining 
higher-order correlation functions of a spin- 2 field such as 
the weak lensing shear or the CMB polarization by using 
the "center of mass" of the configuration as the origin from 
which the components of the shear are defined. In princi- 
ple for an iV-point function there are 2^ different statistics 
of the shear. We have shown that these statistics can be 
divided according to their behavior under parity transfor- 
mation and that one does not expect a cosmological signal 
in the odd ones. 

In order to gain intuition about the behavior of these 
statistics we calculated the four three-point functions un- 
der some simple assumptions. We calculated analytically 
what would be expected to be produced by an ensemble 
of singular isothermal spheres. We showed that C^~^~^ is 
positive and expected to carry the bulk of the signal for 
triangles that are not too elongated. For elongated trian- 
gles we showed all four three point functions have similar 
values but some of them are positive and others are neg- 
ative. If two of the points are very close to each another 
and C^"*"^ carry most of the signal. 

We estimated the three-point functions in the context 
of the halo model using the contributions from the one- 
halo term, which should be a reasonable approximation at 
small angular scales. We showed that the configuration 
dependence in this case is almost identical to that found 
for the SIS. We estimated that at scales of order of one ar- 
cminute roughly N ~ 10^ triangles are needed to detect a 
signal above the noise produced by the intrinsic ellipticity 
of the galaxies. 

Clearly both the fact that we restricted ourselves to 
the one-halo term and that we assumed the halos to be 
spherical will affect the behavior of the three-point func- 
tion. A detailed study using numerical simulations will 
be needed to improve upon the calculation presented here 
(Benabed, Scoccimarro & Zaldarriaga, in preparation). 
The potential rewards of detecting a non-Gaussian signal 
in the shear maps are enourmous and the tantalizing de- 
tections reported so far (Bernardeau, MeUier & Van Waer- 
beke 2002b) make this a very exciting time to study these 
issues in detail. 

When this paper was under completion, a similar pro- 
posal for calculating the shear three-point function was 
put forward by Schneider & Lombardi (2002). 
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